Tunable quantum dots in bilayer graphene 
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We demonstrate theoretically that quantum dots in bilayers of graphene can be realized. A 
position-dependent doping breaks the equivalence between the upper and lower layer and lifts the 
degeneracy of the positive and negative momentum states of the dot. Numerical results show 
the simultaneous presence of electron and hole confined states for certain doping profiles and a 
remarkable angular momentum dependence of the quantum dot spectrum which is in sharp contrast 
with that for conventional semiconductor quantum dots. We predict that the optical spectrum will 
consist of a series of non-equidistant peaks. 
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Two-dimensional (2D) carbon crystals, such as single- 
layer and bilayer graphene, have been the subject of in- 
creasing interest due to the unusual mechanical and elec- 
tronic properties, which may lead to their use in novel 
nanoelectronic devices. The relativistic-like properties of 
carriers in single-layer graphenea>2a2a£ai result from the 
gapless and approximately linear electron spectrum near 
the Fermi energy at two inequivalent points of the Bril- 
louin zone. The charge carriers in these structures are 
described as massless relativistic fermions and are gov- 
erned by the Dirac equation. In contrast, for a symmetric 
graphene bilayer the spectrum is parabolic at the vicinity 
of the K points. 

Among the unusual properties of single-layer graphene 
is the perfect forward transmission across potential bar- 
riers, known as Klein tunneling^, which is related to 
the absence of a gap in the carrier spectrum. This ef- 
fect prevents the electrostatic confinement of charged 
particles and thus the realization of quantum dots. Re- 
cently, alternative strategies have been proposed to con- 
fine charged particles by using thin single-layer graphene 
strips^ or non-uniform magnetic fields^. Here we pro- 
pose a novel approach by considering a bilayer graphene, 
in which a charge imbalance between the layers gives rise 
to a gap in the spectrum, that can be used to create 
potential barrier a 11 ! 12 . Such bilayers of graphene can be 
obtained, e.g., from a graphite crystal by micromechan- 
ical cleavage^. A recent report described the synthe- 
sis of bilayer graphene sheets by graphitization of silicon 
carbide (SiC) surfaces in which the equivalence of the 
two graphene layers is broken by their interaction with 
the SiC substrate as well as by doping one of them with 
potassium atoms«. 

These recent experimental progresses raise the possi- 
bility of introducing a position-dependent modification of 
the spectrum at the Dirac point by changing the potas- 
sium density at different regions of the bilayer graphene 
sheet or by using microstructured gates. In this letter 
we propose a position-dependent potassium doping to 
manipulate the band structure of bilayer graphene in or- 
der to create nanometer-scale quantum structures such as 
quantum dots. Semiconductor quantum dots have been 



intensively investigated both theoretically and experi- 
mentally. Up to now, only quantum dots on graphite 15 
have been obtained, with sample thicknesses of several 
nanometers, i.e., of the order of 10 stacked graphene lay- 
ers. Theoretical calculations show that the properties of 
a bilayer graphene can be significantly distinct from such 
multilayer systems^. In the present work we demon- 
strate that the electronic states of such quantum dots can 
differ significantly from the usual semiconductor-based 
dots. 

The crystal structure of an undoped bilayer of 
graphene is that of two honeycomb sheets of covalent- 
bond carbon atoms coupled by weak Van der Waals 
forces. To each carbon atom corresponds a valence elec- 
tron and the structure can be described in terms of four 
sublattices, labelled A, B (upper layer) and A', B' (lower 
layer). Considering only nearest-neighbor hoping, the 
Hamiltonian of the system in the vicinity of the K point 
is given, in the continuum approximation, byii 



H = H Q + (AU/2)t z 
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where t w 400 mcV is the interlayer coupling term, ir = 
VF(j>x + iPy), P = iPxiPy) is the 2D momentum operator, 
v F w 1 x 10 6 m/s, U = (Ui + U 2 )/2, AU = U\ — U 2 
and U\, U 2 are the potentials at the two layers, which 
reflect the influence of doping on one of them and/or the 
interaction with an external electric field, i.e., the gating. 
The operator r z assigns a positive (negative) sign to the 
upper (lower) layer labels and is defined as 
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with I denoting the 2x2 identity matrix. The eigen- 
states of Eq. (1) are four-component spinors \& = 



FIG. 1: Lower panel: schematics of the quantum dot layout, 
showing the doped graphene layer (brown) and the quantum 
dot region (light brown). Upper panel: the induced parabolic 
potential and the energy levels of the confined states. 



[ipA , tpB , ipB' , V"A'] T , where ipA,B (i>A',B') are the enve- 
lope functions associated with the probability amplitudes 
at the respective sublattice sites of the upper (lower) 
graphene sheet. One can also define an inversion op- 
erator X as 



X = 



I 

1 



(4) 



which acts on \P by changing A to B' and B to A'. It 
is seen that for finite AU, the Hamiltonian (1) does not 
commute with X; this results from the lack of inversion 
symmetry in the system. 

For constant potentials U% and U 2 , the single-particle 
spectrum consists of four bands that correspond to the 
solutions of the equation [s F — (E— Ui) 2 ][s 2 F -(E-U 2 ) 2 } — 
t 2 (E - U X ){E - Z7 2 ) = with s F = hv F k, given by 



E+(k) = Uo + (l/2)[(t±r) 2 + n] 1/2 , 



E ± (k) = [/ -(l/2)[(i±r) 2 + fi] 



1/2 



(5) 
(6) 



where T = [t 2 + As 2 , + A(s 2 F /t 2 )AU 2 ] 1/2 and fi = 
[l - 4s 2 F /t 2 ]AU 2 . Note that for k = the spectrum 
shows a gap at k = of size E±(0) - EZ{0) = \ AU\ and 
the system becomes a narrow-gap semiconductor. Figure 
2 shows the spectra for XJ\ = U 2 = (dashed red curve), 
Ui = 200 mcV and U 2 = (solid black), and U x = 
—U 2 = 100 meV (dot-dashed blue curve). The existence 
of such a gap allows the creation of nanostructures based 
on a position-dependent doping, which in turn can cre- 
ate a position-dependent potential difference between the 
upper and lower layers of the structure^. Different po- 
tentials Ui and U 2 can be created on defect-free graphene 
by controlling the deposition of potassium atoms on the 
upper layer or by applying gate electrodes^. 

The next step is to use a non-uniform doping, or mi- 
crostructured gate electrodes, that will lead to a position- 
dependent potential difference between the layers. In the 
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FIG. 2: Spectrum of a graphene bilayer with Ui = U2 = 
(dashed curve), Ui = 200 meV and U2 = (solid curve), and 
Ui — ~U2 = 100 meV (dot-dashed curve). 



experimental system studied in Ref. 14 there was a po- 
tential difference AU between the graphene layers even 
in the absence of doping, due to their interaction with the 
SiC substrate. In that case the potentials were equalized 
when the number of doping electrons per unit cell was 
70 = 0.0125. For the sake of simplicity we assume a con- 
stant interlayer hopping term t though in reality t may 
depend on the doping level. In the present case a small 
change of the value of t will not affect significantly the 
results for low energies. The experimental results show 
that AU scales approximately linearly with doping over 
a wide range of potassium concentration levels. 

A position-dependent, circular-symmetric profile 
7(1*) = 7(/o) gives rise to a 2D circular-symmetric 
quantum dot of radius R on a bilayer graphene, shown 
schematically in Fig. 1, with the doping level ranging 
from 70 at p = to a maximum value 7m for p > R. In 
that case the four-component spinors ^ are given by^ 



( Mpy m6 \ 

i(t)B(p)e l{m - 1)e 
4>B>{p)e me 



(7) 



V*Mp)e <(m+1)fl / 
where m is an integer. ^ are eigenstates of the operator 



J *- L ' + ^ T "-^\0 -a, 



(8) 



with eigenvalue to; L z is the angular momentum opera- 
tor and o~ z the Pauli matrix. This represents the total 
angular momentum of the carrier and the layer index 
operator, which is associated with the behavior of the 
system under inversion. 

The potential difference between the layers is AU = 
U(p) and subsequently the equations for the radial part 
of the spinor functions are 



[d/dC + m/(\(j) A = Sltfa, 



(9) 



FIG. 3: Lowest energy states of a graphene bilayer quantum 
dot for m = (blue solid line) , m = 1 (black dotted line) and 
m = — 1 (red dashed), for (a) Um = 50 meV and U2 = and 
(b) Ui = -U a and [7m = 100 meV. 



[d/d( - (m - 1)/C]0s = -fli^i + (10) 
[d/dC+(m + l)/C}(t>A' = -^(t>B'+t'(j)A, (11) 



[d/c!C-m/C]0 



B' 
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where C = p/R, ^1,2 = e — "1,2, "1,2 = Ui^R/hvF, 
e = ER/hvF, and t' — tR/hvF- The mil that ap- 
pear in Eqs. (10) and (11) are a direct consequence of 
the presence of a geometric phase in the system (i.e. a 
full rotation of the system introduces a phase difference 
of 2ir between the sublattice sites A' and B). In addi- 
tion, the large magnitude of the layer coupling term t', 
in comparison with the other energy terms, implies that 
the 4>b and 4>A' must have significantly higher amplitudes 
than <pA and <pB'- These two facts have important con- 
sequences for the energy spectrum and the shape of the 
resulting probability density functions, as we show below. 

For the numerical calculations we took a truncated 
parabolic profile AJ7 = UmC 2 > for £ < 1 and AU = Um 
otherwise. The energies of the lowest energy bound states 
are shown in Fig. 3 as function of the radius of the dot. 
Figure 3(a) shows the case U2 = 0, Um = 50 meV and 
Fig. 3(b) the case Ux = -U 2 , U M = 100 meV. In both 
cases the spectrum is proportional to the inverse of the 
dot radius and each set of energy levels are approximately 
equally spaced for fixed m. The results correspond to the 
lowest energy levels for m = (blue solid lines), m = — 1 
(red dashed), and m = +1 (black dotted lines). An im- 
portant difference between (a) and (b) is the absence of 
hole bound states in (a), since there is an overlap between 
the hole band states in the well and the corresponding 
band in the barriers, whereas in (b) both electron and 
hole confined states are present. We found the remark- 
able result that the ground state is almost degenerate for 
m = and m = il in both cases (a) and (b). This 
turns out to be a consequence of the geometric phase 
which shifts the angular momentum labels in Eqs. (10) 
and (11). In comparison to the results for semiconductor 
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FIG. 4: Panels (a)-(c): Radial dependence of the spinor func- 
tions for a parabolic quantum dot on a doped bilayer graphene 
for m — 0, and m = il. Panel (d): probability density for 
each case. The parameters are R — 50 nm, Ui — 50 meV and 
U 2 =0. 




FIG. 5: Energy spectrum as a function of angular momentum 
label m for R = 50 nm, and Ui = —U2 = 50 meV. 



quantum dots we find that in a graphene bilayer there 
is a peculiar energy shift between states with positive 
and negative angular momentum even in the absence of 
an external magnetic field. This is a consequence of the 
lack of inversion symmetry in the system when U± Ui- 
The shift is larger for the low-lying states, with the +m 
and — m states becoming nearly degenerate as the en- 
ergy increases. A similar, albeit weaker, shift is found 
for |m| > 1. 

The radial part of the spinor functions for some lowest 
energy states is shown in Fig. 4 for to = (a), m = — 1 
(b), and m = 1 (c). Panel (d) shows the probability den- 
sity for Tfi — (solid black curve), to = —1 (blue 
dot-dashed curve) and m — 1 (red dashed curve). All 
results were obtained with U± = 50 meV, R — 50 nm, 
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FIG. 6: Oscillator strengths for transitions between the lowest 
energy levels (m, n) in quantum dot of radius R — 50 nm with 
Ui = -U 2 = 50 meV. 



and U2 = 0. In contrast to a semiconductor quantum 
dot, the present results show that the probability den- 
sity for to = has a maximum at p ^ 0. This is due 
to the fact that, for to = 0, the high-amplitude spinor 
components (i.e. <pB and 4>A') vanish at the origin. In 
contrast, for m = 1, which is here the ground state, one of 
the high-amplitude spinor functions is allowed to assume 
non-zero values at the origin. In other words, for to = 1, 
an electron encircling the dot adds an angular momentum 
phase of 2ir which can be cancelled by a Berry's phase 
of — 2ir. The shift between the low- lying energy levels of 
positive and negative momentum labels is evident in Fig. 
5, which shows results for a quantum dot with U2 = —U\ 
and R = 50 nm. The results is this case show that the en- 
ergy eigenvalues are related by E(m) = —E(—m). This is 
consistent with Eqs. (9)-(12) and in sharp contrast with 



the relation E(m) = E(—m) satisfied in a usual quantum 
dot with parabolic confinement. A similar behavior was 
observed for other potential profiles, such as the square 
quantum dot and the tangent hyperbolic potential. 

Figure 6 shows the oscillator strengths for the allowed 
transitions (to, n) — > (to ± 1, n') between the three low- 
est energy states (i.e. to = and ±1) within the energy 
window of < AE < 9 meV, in a quantum dot of ra- 
dius R — 50 nm with U\ = — U2 = 50 meV. Since the 
energy levels for different to are not equally spaced, the 
results show the matrix elements at several values of AE. 
This differs considerably from the case of the parabolic 
quantum dots in semiconductors, in which the transitions 
occur at a single energy. These transitions can be probed 
by means of far infrared spectroscopy measurements. 

In summary, we proposed a novel approach to electri- 
cally confine charge carriers in 2D quantum dots in bi- 
layer graphene. The quantum dots are created through a 
position-dependent doping which breaks the equivalence 
between the upper and lower layers. A further tuning of 
the quantum dot can be realized through the use of gates. 
The numerical results show that the degeneracy of the 
positive and negative momentum states of the dot is lifted 
even in the absence of an external magnetic field. This 
result differs fundamentally from that for conventional 
semiconductor quantum dots and arises due to the lack 
of inversion symmetry caused by the doping. The sys- 
tem can be realized experimentally by a suitable choice 
of doping levels or with the application of electric fields. 
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